Two-Body Bound States in Light-Front Dynamics 



J.Carbonell (a) , M. Mangin-Brinet (a) [|, V.A. Karmanov (b) 

( a ) Institut des Sciences Nucleaires, 53 avenue des Martyrs, 38026 Grenoble Cedex, France 

( b ) Lebedev Physical Institute, Leninsky Prospect 53, 119991 Moscow, Russia 



O 
O 

<N 

two scalar particles in the Wick-Cutkosky model and of two fermions interacting via the 
pL( | usual OBEP ladder kernels. 
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Abstract 

We present the main features of the explicitly covariant Light-Front Dynamics formal- 
ism and a summary of our recent works on this topic. They concern the bound states of 



<N : 1 Motivation 

CSj ■ The recent measurements performed at CEBAF/TJNAF on the deuteron structure functions 
and tensor polarization up to momentum transfer values yet never reached, was the starting 
motivation for a series of works in relativistic dynamics. Indeed the measurement of the elec- 
tromagnetic form factors at values Q 2 ~ 6 (GeV/c) 2 , much greater than the nucleon mass, 
lets no hope for any attempt to describe this physics without relativity. There are even very 
few chances that meson-nucleon dynamics could still be used while the distances involved are 
smaller than the nucleon itself Rn ~0.86 fm. A recent and complete reviews of both the 
experimental and theoretical works can be found in [jl], |2| and references therein. 

This motivation has always been present in theoretical physics since the first days of 
quantum mechanics. It relies on the need for having a proper relativistic wavefunction of 
simple systems - an essential ingredient in the NN, qq, qqq, . . . physics - and for clarifying the 
kind of modifications that Lorentz invariance brings into the non relativistic dynamics. 

Our first contribution in this field consisted in a series of papers [§, f§, fj], |§ which, in 
the framework of the Explicitly Covariant Light-Front Dynamics (ECLFD), led us to predict 
the deuteron structure functions and polarization observables. These calculations, though per- 
formed in a perturbative way over the wavefunctions of the Bonn model 0, turned to be quite 
successful in the deuteron case and have recently found a natural application in describing the 
NN correlation functions in nuclei (see Prof. Antonov contribution to this School and ||). 
The success encountered in these first works and some clear advantages of the ECLFD itself 
in describing relativistic composite systems, convinced us to go deeper inside this approach 
both in developing formal aspects of the theory and in studying in detail the properties of its 
solutions. This effort has up to now resulted in studying the bound states solutions of two 
scalar and two fermions systems in the ladder approximation. Doing so, we have in mind the 
description of "genuine relativistic systems" (those for which B ~ m, (3 ~ 1) as well as the a 
priori non relativistic ones when probed in relativistic regions (e.g. deuteron). 

1 Now at Geneva University and CERN, Switzerland 
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In this contribution, we first present the leading ideas of the ECLFD formalism (Section 
HD and then summarize the main results of our recent works f9|l- ||18| on this topic, both in the 
scalar (Section |3|) and fermionic (Section |]) case. We put special emphasis in the following 
items: i) the size and nature of relativistic effects ii) the comparison with the non relativistic 
solutions and other relativistic approaches iii) the problem of constructing the non zero angular 
momentum states and iv) the solution of the OBEP models for fermions and their stability 
properties with respect the cutoff. 



2 The formalism 

Light-Front Dynamics is an hamiltonian formulation of the quantum field theory, specially well 
adapted to the description of relativistic composite systems. In its explicitly covariant version, 
proposed by V.A. Karmanov [19| and recently reviewed in || |20| , the state vector is defined 



on a space-time hyperplane whose equation is given by uo ■ x = a with uo 2 = 0. Wavefunctions 
defined as the Fock components of the state vector - are the usual formal objects of this theory 
and are directly comparable to their non relativistic counterparts. 

As in any hamiltonian formalism, the dynamical equations are provided by the eigenvalue 
equations of the two Casimir operators of the Poincare group, respectively the square of the 
four-momentum P M and of the Pauli-Lubanski four-vectors 

P 2 \m> = M 2 \m> (1) 
£ 2 |^> = -M 2 J(J+1)\^> (2) 

These operators are quadratic forms of the ten generators j-P^, J p(T | of the Poincare group: 

p2 = 

Different ways of constructing the generators set would lead to different relativistic theo- 
ries 0] which will be furthermore strictly equivalent if they were solved in their full complex- 
ity. This choice is provided by integrating the Noether conserved currents over a suitable 3D 
space-time surface. We remind that Noether theorem associates to each continuous space-time 
transformation 

x — > x = Ax + a (3) 

letting invariant a lagrangian density 

£{x) = £ {x) + £ int {x) (4) 

a 4-current satisfying the continuity equation 

dpf = (5) 

This latter is calculated from a particular £(x) according to very precise recipes. Integrating 
over an arbitrary 4D volume V (figure H), equation (^) is equivalent to assert the nullity of flux 
across its boundary dV 

' fdS, = (6) 
av 
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Figure 1: Arbitrary volume V for integrating continuity equation 



where dS^ denotes the normal to dV. 

Equation @ leads to a conservation law only if V is such that "lateral flux" vanishes. In 
the usual approach of dynamics - denoted "instant-form" in Dirac classification [|2TJ - V is the 
space-time volume limited by the hyperplane t=cte (see figure §). In Light-Front Dynamics - 
denoted "front-form" in |2l| - V is limited by two space-time planes u ■ x = a (see figure |||). 
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Figure 2: Flux in instant-form of dynamics Figure 3: Flux in front-form of dynamics 



Let us illustrate the procedure in case transformation (|3J) is a translation 

/ , 

Xy ► Xy Xy "T Ojy 

Noether theorem provides a set of four 4-currents {T^ v } u=0 3 given by 

dC 



rpflU 



d v <t> - g^C 



and satisfying 

d^(x) = 

They constitute the energy-momentum tensor. Integrating the continuity equation over an 
arbitrary space-time V one gets 

J d 4 x d u T" v = 

and choosing V like in figure ||| gives 



f dxT 0u (x) 


+ 


J dxT 0u (x) 


= 




upper 




lower 



That means the four quantities 



£"(to) 



dxT 0v (x 



are independent of the hyperplane t — to position. 

Alternatively, by choosing space-time planes ui ■ x = a, one obtains the cu-dependent 
generators 

P"(tr) = f T^{x)5{u -x-a) Lu u d A x (7) 



with 

d a P"(<T) = 

Equation (0) is the starting point of the Light-Front Dynamics. 

In the instant-form approach, the spatial components depend only on the non interacting 
part of the lagrangian and are called kinematical. The zero-th component is the interaction 
Hamiltonian and governs the dynamics 



P 



Hint 

Po 



In LFD, all generators split into a free part (Pq) and an interaction-dependent one (Pf^t) 



P»(a) = PX + PLM 



with the interaction dependent part being along the tu-direction 



Ptnt{ a ) = ^ / H int (x)5(uj ■ x - a)d A x 



This can be turn to profit to obtain a maximum number of seven non interacting generators, 



an advantage that was already noticed by Dirac in his seminal paper 0] about the different 
forms of relativistic dynamics. 

Once obtained the generators, (0) provides the dynamical equation determining the mass 
of the system M 2 . After some algebra, one is led to 



(M 2 - P 2 ) 


* >= 2P • w j 


' Hi n t{oJT) exp(-zcrr) dr 


* > 



(8) 
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where Tli n t(k) denotes the Fourier transform of the hamiltonian density 



H in t{k) = J Hi n t(x) exp(ik-x)d 4 x 



We will consider hereafter the a-stationary solutions and set a = in 
states towards a is given in |T8| . 

The state vector is decomposed into its Fock components ^ a /3 f\ 



The evolution of 



* > = f • • • d 4 k a d 4 qi . . . d% 

a/3 J 

Vupih, . . . , k a , gi . . . q p ) 4 X . . . 4 q 6j 1 . . . b\ | > 



(9) 



which are the wavef unctions. We can consider the set ^n af} = {^a/?} as the components of 
an infinite dimensional vector ^ = (^i, ^2, ^3, ■ ■ ■), coupled to each other via the interaction 
operator / 7i irii (u;r) dr. Equation (||) is thus an infinite system of coupled channels. 



(M 2 
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If we restrict ourselves to the two- = {^20}) an d three-body (^3 = {^21}) wavef unctions, 
we obtain a system of two coupled equations for and \&3 which constitutes the ladder 
approximation. By expressing ^3 in terms of ^2, one gets an integral equation for ^2 with 
energy dependent kernel. 

The construction of non zero angular momentum states in LFD is made difficult by the 
fact that generators of spatial rotations contain the interaction. Consequently the eigenvalue 
equation for S 2 and S3 are non trivial dynamical equations with a level of complexity higher 
than the mass equation (|S]) one. To circumvent this difficulty we have made use of a kinematical 
angular momentum operator 

J = —i \ k x — + n x — 
\ dk dnj 

whose eigenstates can be constructed by using the standard angular momentum algebra. The 
physical solutions are then obtained as a linear combination of the kinematical operator A 2 = 
(J • n) 2 eigenstates, satisfying the so called "angular condition". This method was shown to 
restore - at least in a simple model - the rotational invariance of the theory to a high degree 
of accuracy [17|. Its explanation requires long technical developments and will not be detailed 



here. The interested reader could find the first published results in |jlQ|, O, |T7 



An important property of Light-Front Dynamics is the fact that the vacuum of the theory 
is decoupled from any other Fock state despite the interaction. This "emptiness of the vac- 
uum" results into the absence of vacuum fluctuations diagrams and considerably simplifies the 
evaluation of certain processes. 



2 We consider two interacting fields with corresponding creation operators aj, and 
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An appreciable advantage of this formalism with respect to other relativistic approaches, 
is the clear link with the usual non relativistic dynamics: LFD wavefunctions have the same 
physical meaning of probability amplitudes than their non relativistic counterparts. 

Among the drawbacks we note - apart from the psychological barrier of using a non 
conventional formalism - the appearance of contact (instantaneous) interactions for non scalar 
constituents and two direct consequences of solving the dynamical equations in a truncated 
Fock space: the uo dependence of the on-shell approximate amplitudes and the fact that the 
different projections of the non-zero angular momentum states along the u direction are non 
degenerate (the so called violation of the rotational invariance). 

It could be of some help to say few words about the words in what concerns the Light- 
Front world. Historically it appeared first in Dirac classification with the particular choice 
of front surface t — z — 0. It reappeared in Quantum Field Theory in 1966, under the label 
"Infinite Momentum Frame" calculations, when it was realized J23| that boosting the frame 
in which calculations of invariant amplitudes are made, gave easier results. In particular all 
diagrams related to vacuum fluctuations automatically vanished. Later on, it was shown p4 



that this procedure was equivalent to use the LFD variables z ± t in the theory from the 
very beginning and constitutes now the so called standard LFD, developped by many authors 

n, m, m, m s m m in n- 



The difference with respect to our approach is in the words " Explicitly" and " Covariant" 
which we will try to justify. To choose a front surface in the form to ■ x = o is not only a 
mathematical " delicatesse" but a way to carry everywhere in the theory the w-dependence 
of the formal objects in an explicit way. It has several advantages, all related to the fact 
that to is a four vector with well defined transformation properties. For instance it provides 
explicitly covariant expressions for the on shell amplitudes, a property which is often hidden in 
the standard LFD formulation. The latter one is recovered by fixing the value uj = (1, 0, 0, — 1) 
but this value is associated to a particular reference frame and it is not valid in any other one, 
what implies a non covariant formulation. Another advantage of the covariant formulation is 
the fact that, because the cu-dependence is explicit, it can be controlled at will and removed 
when calculating observables (e.g. in the form factors f35fl). 

If by covariance of a theory, we understand the fact that is based on a set of generators 
satisfying the Poincare algebra, then the standard and the explicitly covariant approaches would 
be both covariant if they were solved in their full complexity. This is never the case in practice, 
in particular due to the Fock space truncation, and both formulations can lose this property 
when approximate schemes are used. 



3 Two Scalar particles 

It is always instructive to start by considering a simple model which could provide us with 
some physical insight with a lower formal cost. We have thus considered the Wick-Cutkosky 
(WC) model |3B| describing two scalar particles (m) interacting by a scalar exchange (/i) with 
a Lagrangian density Ci = g(p 2 X- We are interested in bound states of two particles with 
momenta k\ and /c2- In the reference frame k\ + k2 = the ECLFD equation reads 

— * in ^ r d^h'^ —* —> — * 

[A{k 2 + m 2 ) - M 2 }^(k, ti) = -— —V(k, k>, h, M 2 )V{k>, h) (10) 

27P 1 J €k> 
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where k—ki, e& = \fk 2 + m 2 , co> = (l,n) in this particular frame with v? 
mass of the system, related to its binding energy B by M = 2m — B. 



1 - and M the total 



The interaction kernel is[] 



V(k, 9, h, M 2 ) 



4ira 



Q 2 + ^ 



with 



Q 2 = (k - k'f - (n-k)(n-k') 



(efc' — £k) x 



+ \ei + et- 



M 2 



n-k' n-k 



Relativistic effects are all included in the second and third n-dependent terms of Q 2 . By 
formally setting n = in this expression - though n 2 = 1 ! - one recovers the non relativistic 
Coulomb kernel. By furthermore assuming B << m, the kinematical left hand side term of (|10|) 
becomes 4m(k 2 + B) and the ECLFD equation coincides with the non relativistic Schrodinger 
equation in momentum space for Coulomb interaction^. 

Equation fllOD was solved ||, [18|] for several values of J and ji. For S- waves, the solution 
\I/ is a scalar function depending on two scalar arguments ^ = ^(k,n ■ k), i.e. the LFD S-wave 
wavefunction has an angular dependence! For non-zero angular momentum states the solution 
was obtained on a basis of eigenfunctions of the A 2 operator, discussed in the previous section. 
In a truncated Fock space, these eigenfunctions were found to have different M 2 values [H| [12 



- a result also noticed in [57] - and a method was proposed to restore in the two-body sector the 
rotational invariance. Our main results concerning S-waves are summarized in what follows. 
The numerical values given hereafter all correspond toh = c = m = l units. 



Size and nature of relativistic effects 



The size and nature of relativistic effects in the binding energies can be seen in figure |]. 
The LFD results for fi = are compared to the non relativistic ones and to those given by the 
Klein-Gordon equation. One can see a quick departure from Schrodinger results and a clear 
repulsive character - smaller binding energy - contrary to Klein-Gordon (and Dirac) equations 
which produces attractive corrections to all orders. LFD energies are also compared to a first 
order perturbative calculations - equally valid for Bethe-Salpeter (BS) equation - which have 
the form @ 

ma 2 ( 4 , 1 \ . . 

B pert = — 7- l--cdog- 12 
4 V 7r aj 

and turn to be relevant until a ~ 0.3. 

It is worth noticing that sizeable relativistic effects are observed in a system for which 
both the binding energy and the average momentum are small. At a = 0.3, for instance, they 
account for a 100% effect in the binding energy whereas <fc 2 > w 2%. 

We would also like to notice that the large differences between the ECLFD and the 
non relativistic solutions are not of kinematical origin. In order to disentangle the different 
contributions to the relativistic energies B, equation (|10D has been formally written K^> = ^V^/ 

3 Strictly speaking, the Wick-Cutkosky model corresponds to the /i = case. We consider here a trivial 
extension 

4 We note a missprint in page 465 of pdj| where the expression Amk 2 + 2mB was instead written 
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Figure 4: LFD energies (solid) for \x = as a Figure 5: Contribution of different relativistic 
function of a, compared with non relativistic, corrections (see text for details) 
Klein-Gordon and perturbative results 

and we have considered several approximations of it. The results are displayed in figure We 
have first considered the case of a non relativistic kernel V - i.e. put n = in ( |TT|) and e = m 
- together with in curve a the non relativistic kinematics K = 4m(k 2 + B) and in curve a' the 
relativistic one K = 4(k 2 + m 2 ) — M 2 . Curves b and b' are obtained in the same manner but 
putting e = Vk 2 + m 2 . The last one corresponds to the full ECLFD equation. These results 
show that - at least in what scalars are concerned - the main effect of a relativistic formalism 
is not in the kinetic energy but in the interaction kernel. The kinematical term K has a very 
little influence on B, which furthermore goes on the opposite direction, whereas e and V kernel 
contribution are both essential. One can thus conclude that kinematical corrections alone are 
not representative of relativistic effects. Even when they are included in the kernel, e.g. through 
efc, they can only account for half the effect. 

Comparison with Bethe-Salpeter equation 

The comparison with Bethe-Salpeter equation is done in figures ^| and [7|. Figure |^ rep- 
resents the LFD B(a) curves (solid line) for different values of /i. They are compared with 
those provided by BS equation (dashed-line) in the same ladder approximation, whose kernel 
incorporates higher order intermediate states. Their results are seen to be close to each other. 
This fact is far from being obvious - specially for large values of coupling constant - due to 
the differences in their ladder kernel. A quantitative estimation of their spread can be given by 
looking into an horizontal cut of figure || i.e. calculating the relative difference in the coupling 
constant (ojlfd — c^bs) I '®lfd for a fixed value of the binding energy. The results, displayed in 
figure |7| for B = 1.0, 0.1, 0.01, show that relative differences (i) are decreasing functions of \x for 
all values of B (ii) increase with B but are limited to 10% for the strong binding case B = m 
which involves values of a > 5. This indicates the relatively weak importance of including 
higher Fock components in the ladder kernel even for strong couplings, as was noticed in . 
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Figure 6: B(a) for different values of [i in LFD Figure 7: Differences in the coupling constant 
(solid) and BS (dashed) approaches as a function of \i for fixed values B. 



Weak binding limit and non relativistic solutions 

An interesting result was obtained by studying the weak binding limit of B(a). We found 
that, except for the case of zero mass exchange, relativistic and non relativistic solutions differ 
even when describing zero binding energy systems. These results are displayed in figures |8] 
and H BS solutions, also included, display the same behaviour, indicating that this is not 
a pathology of the Light-Front formalism but seems rather a general feature of consistent 
relativistic theories. A system bound by a massive exchange would be described by different 
parameters whether one uses LFD (and BS) or Schrodinger approaches, no matter how small 
the binding energy will be. We concluded from that to the non existence of non relativistic 
limit for /i ^ 0. 

We would like to mention here that an important part of these differences were taken into 



account by including in the Schrodinger equation energy dependent interactions |40 . 



The scalar deuteron 

A straightforward application to this model, and a way to evaluate the modifications 
induced by a relativistic treatment of a semi-realistic model, was done by building a scalar 
model for deuteron. By adding a "repulsive scalar" term[| one gets a relativistic version of the 
Nucleon-Nucleon MT potential []39| . 



v = v R + v A 



X- 



(13) 



In the non relativistic limit, A#, /iR, A^, ha are adjusted to reproduce a realistic deuteron wave- 
function with B = 2.23 MeV and acceptable NN scattering parameters. If one inserts kernel (|P3|) 

5 One should note that such a repulsive interaction cannot arise from a scalar exchange but can mimic some 
fcrmionic interactions 
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Figure 8: Zero binding energy limit of LFD Figure 9: The same study for different val- 
and BS equations (solid) compared with non ues of \i. LFD and BS curves are not distin- 
relativistic solutions (dot-dashed) for \i = 1 guished 

in the LFD equations, the binding energy is shifted from B = 2.23 to B = 0.96 MeV, a dramatic 
repulsive effect. One can recover the physical value for B by shifting A# = 7.29 — > Xr = 6.60 - 
all other parameters being unchanged - what makes a 10% difference in a coupling constant. 
The spectacular change in B is partially due to the small value of the deuteron binding energy. 
It illustrates however well the difficulty to unambiguously determine a coupling constant in the 
strong interaction physics, even for systems widely considered as being a priori non relativistic. 

Wavef unctions 

The main modification when comparing LFD and NR wavefunctions with the same cou- 
pling constant is due to the change in the corresponding binding energies. A first requirement 
to compare them is thus to re-adjust one of the a values (we choose the NR one) in order to 
deal with states of equal B. Figure [10] shows the comparison of these wavefunctions, squared 
and integrated over the angular dependence n • k = cos#, for ji = and a moderate coupling 
a = 0.5. One then has 5^^=0.0267 and -Bjvj?=0.0625 which has been re-scaled to the LFD 
value with a = 0.327. One can see small deviations (~5%) at k — 0. In the relativistic region, 
the LFD solutions are systematically smaller than the NR ones and their differences increase 
reaching a factor 3 for k = 2, despite the moderate values of B and a. When dealing with highly 
relativistic systems - like e.g. mesons in the constituent quark model - both descriptions differ 



even for small values of k. Figure 11 shows the case Bifd=0.S4: and a = 5 where the squared 



wavefunctions at k = are ^70% different. 
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Figure 10: LFD and Schro dinger wavefunc- Figure 11: Similar comparaison in a highly 
tions, squared and integrated over 6 for /i = 0. relativistic case (a = 5) 

4 Two fermions system 



We have also obtained |L3[ [T^] recently the LFD bound state solutions for a system of two 
fermions interacting with the usual - scalar (S), pseudoscalar (PS), pseudovector (PV) and 
vector (V) - OBEP Lagrangians: 



C 



9pv 



C pv = ^ lbl d^ 
2m 

As in the scalar case, calculations were performed in the ladder approximation. The interaction 
kernel is provided by the Born scattering amplitude M. and - according to the Light-Front graph 
technique - has two contributions which differ from each other by the light-front time ordering. 
Its analytical structure has the form 

M = {A 1 (k' 1 ,ki,r')+A 2 (k 1 ,k' 1 ,r)} [u^k'^u^h)] [ua! 2 {k 2 )0 2 u a2 {k 2 ) 

where Oi are vertex operators depending on the type of coupling (Oj = 1, 275, . . .) and Ai, A 2 
scalar functions. In "center of mass" variables, i.e. boosting to a reference frame in which 
k\ + k 2 = 0, the scalar part can be written as 



Ax(k,k',n, M)+A 2 (k,k',h,M) 



with Q given by (|lT|) and the total amplitude reads 



Q 2 + V 2 



Q 2 + ti 2 



u^k'^OxUa^kx) u(k' 2 ) a < 2 d 2 u(k 2 ) 
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M 




+ 



q t 



u(k'1) 



u(k1) 



k'1 +cot" = k1 + q -hot ' 




u(k'1) u(k1) 
k'1 + q = k1 +cot -km; 



As in the scalar case, the two-fermion LFD wavefunctions $ are Fock components of state 
vector (H). Their spin structures depend on the quantum numbers of the state. For J n = + 
they have the form 

$ SlS2 (A;i, k 2 ,p,u) = u S2 (k 2 ) <p lbU c u Sl (ki) (14) 



with 



+ 



2m7 /J w fI 4m 2 



lu ■ p s 

depending on two scalar functions <pi. For J n = 1 + they read 



Ks 2 ( k i, k 2 ,P,u>) = ej(p) u S2 {k 2 ) r U c u Sl {K 



(15) 



in which 



+ 



2m z m uj ■ p 

i ^T5^ fJ.U per 

(ki + k 2 ) u (ki — k 2 ) p uj r7 



(fci - fegVjj 
2m(w • p) 



2m?uj ■ p 



+ 06" 



uj ■ py 



p is the total momentum, Uc = 7270 the charge conjugation operator, e*(p) the polarization 
vector and u = 7^^. They depend on six scalar fonctions fa. 

To get easier links with the non relativistic wavefunctions one uses instead of fa some 
linear combinations denoted /j. For + states, for instance, one simply has 



fi 
h 



2y/2e k 
2y/2k- 



with j\ tending to the usual NR wavefunction and f 2 being an extra component of relativistic 
origin with no counterpart. For 1 + states the combination is more involved || but the fi 
functions have similar properties. Thus f\ tends to the non relativistic S-wave component u$, 
f 2 tends to the D- waves —Ujj and /3_6 are extra components of relativistic origin. 

One can see that in this /-representation, LFD provides a very interesting parametriza- 
tion of the relativistic dynamics. One deals with similar formal objects, depending on similar 
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arguments. In this way, relativity acts by modifying the usual wavefunctions and by introducing 
extra components which becomes negligible - or tends to zero - in the non relativistic regime. 
In $ representation (|T4T - |T5|) the two-body LFD equation reads 



M 2 -{k 1 +k 2 f * 2 ,p,wr) 



m 



E 



2tt 3 ti J 2e fc / 2e 



d 3 k[ d% dT , 5 ^ + k ^_ p _ ur) (16) 



2(u)'p) x (A^, k 2 ,p,ur') x V^a^k'^ k 2 ,p,ujr r ; k h k 2 ,p,ujr) 

Inserting (|1^|15|) in (|16D one is left - after some lengthy algebra - with a set of coupled two- 
dimensional integral equations on the form 



A a (k)(f)a{k, u) =^2 J dk'dv! B a p(k,u,k' \u')4>p{k' ,u') 



(17) 



where u = h ■ k. A a contains kinematical terms and B a p the interaction kernel. For J 71 " = + 
there are two coupled equations whereas for J 71 " = 1 + their number is six but can actually be 



decoupled into 2+4 JT8[ . 

Bij results from integration over tp', the azimutal angle between k and k', of more basic 
kernels k. 



ij 



Bij(k, ii, k', u') 



1 k' 1 r 2 - 



7T" ey Jo 



Kij(k } it, k', u' } if 



, , * F 2 (Q 2 ) df> 



Q 2 + /i 2 2tt 



which have the general structure 

Kij(k, u, k', u, if') = - L a p Tr {Mf"} 



ad 



F is a vertex form factor (FF) depending on Q 2 , L a p is a tensor depending on the kind of 

r/3a 

ij 



coupling, M(! a is a contraction of a product of "7-like" matrices with a tensor IL^ depending 



on the J 71 " of the state. 

= IV S? (k 2 + m) 0% (k' 2 + m) S'] (k'i ~ m) O a x (h - 



ij 



m) 



where k = and ,<S"f are some spin operators |14|, |I5| . 

Kernels Ky are in general calculated numerically. Only for some special cases it is in- 
teresting to deal with analytical expressions, which very soon turned out to be unreasonably 
lengthy. As an example we give below the scalar kernel for J=0 + states: 



K U = 


-a s ii j( 


Kl2 = 


— a s irm 


K 2 i = 


—a s irm 


K-22 = 


-a s 7T {( 



{ (4 - 4) (k sin 9 + k' sin 9' cos tp') } 



2 

where we denote cc. s = |j. Other analytical expressions have also been obtained and can be 



found in [13, 14, M, 18 
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4.1 Results 

In a series of works still in progress we have separately studied, coupling by coupling, 

J 71 = ± , 1 ± states in the loosely bound (B<<m) and ultra relativistic (B ~ m) limits. We will 
restrict ourselves in this contribution to scalar (S) and pseudoscalar (PS) couplings with special 
emphasis in the stability problem - are equations soluble as they are provided by Quantum 
Field Theory or do they rather need some regularization? - and in their comparison with the 
non relativistic solutions. 



2.25 




Figure 12: Cutoff dependence of M 2 in the Fi § ure 13: Cutoff dependence of the coupling 
J = 0+ state, for two fixed coupling constant constant, for J = 1+ states and B = 0.05. 
below and above the critical value. 



4.1.1 The stability problem 

Let us solve the LFD equations in a compact domain [0, k r 



M 2 - A(k 2 



m 



f(k,u) 



dk' 



dv! B(k, u, k', u')f(k', u') 



and ask ourselves what happens with M 2 when k max — > oo. In case the solution exists, we will 
say that equations are stable. In the opposite case, the theory would require some regularization 
procedure, for instance by means of vertex form factors, and it is then pertinent to inquire how 
the solutions depend on them. We have shown that the answer to this question critically 
depends on the type of coupling, on the J n quantum numbers of the state and on the value of 
the coupling constant a. 

Consider first the scalar coupling (Yukawa model) for J=0 + . We found |TJ| the existence 
of a critical coupling constant a c ~ 3.72, below which the solution exists, and above which 



the system "collapses". This behavior is illustrated in figure 12 where we have plotted the M 
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dependence on k max for two different values of a. We showed the asymptotical behavior of 
solutions to be 



with a relation (3{a) provided by an eigenvalue equation that allows a precise determination of 
the critical value /3(a c ) = 0. This property was checked by a direct inspection of the numerical 




Figure 14: Asymptotics of fi for J = + state Figure 15: Asymptotics of fi for J = + state 
with 5=0.05, a=1.096, /i=0.25. The power- with B = 0.5, a=2.48, //=0.25. The power- 
law coefficients are (3\ = 0.82 and /?2 ~ law coefficients are (5\ = 0.55 and /3 2 ~ 



solutions of LFD equations - as it can be seen in figures [14] and [15] - and turns to be very 
accurate. It is worth noticing that - at least in the framework of this model - the coupling 
constant could be "measured" in the tail of the wavefunction! The J=l + states, on the contrary, 



do not present any stable solution without vertex form factors. This can be seen in figure |13 
where the divergency of the coupling constant as a function of k max is displayed for the two 
angular momentum projections of a J=l + state with B = 0.05. One can also remark in this 
figure the non degeneracy of the different projections due to the Fock space truncation; for the 
scalar coupling and moderate values of k max ~ 10, it remains however less than one percent. 

For pseudoscalar coupling, the stability analysis was performed using the same methods 
than for the scalar one fL6| , p~8|| and presents some peculiarities. Quite surprisingly, the equations 
for J=0 + states were found to be stable without any regularization. The results were however 
strange in the sense that they lead to a quasidegeneracy of the coupling constants for binding 
energies which can vary over all the physical range [0, 2m]. One gets for instance a = 49.5 for 
B = 0.001 whereas a = 48.6 for a binding energy 500 times bigger B = 0.5. The origin of 
this anomaly was found to lie in the second channel equation (^22)- ^ nas been understood 
analytically [16|] with a simple model but leads to physically unacceptable results. 

The case of other couplings (PV,V,T) has been examined in [Tj|. Corresponding equations 



results into more singular kernels which all lead to unstable solutions without form factors. 
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4.1.2 Comparison with non relativistic solutions 



We will present the results for J = + states - so with a two-component wavefunction - obtained 
with the scalar and pseudoscalar coupling. In each of them, we will consider a loosely bound 
(B=0.001) and a very relativistic (B=0.5) system. We have fixed for all the calculations an 
exchanged mass \x = 0.15. 
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Figure 16: LFD and non relativistic wave- 
functions for J=0 + state with B = 0.001, 
/i = 0.15 in the Yukawa model 



Figure 17: Comparison between LFD 
(Yukawa model), WC and NR B(a) in J = 
+ state 



Scalar coupling 

For B = 0.001, the LFD coupling constant is aLFD=0.331 whereas the non relativistic 
value is tt7VR=0.323. Like in the Wick-Cutkosky model discussed in the last section - scalar 
particles interacting by scalar exchange - relativistic effects are repulsive. They are responsible 
for only a 3% difference in the coupling constants, whereas in the purely scalar case this 
difference is sensibly greater (avyc=0.364). Wavefunctions - suitably re-scaled - are displayed 
in figure [16]. One sees that component fx dominates over f 2 in all the interesting momentum 
range and that f 2 has a zero at k « 0.25. One can also remark that fx is very close to the 
NR wavefunction in the small momentum region but it sensibly deviates with increasing k; for 
k ~ 1.5 the difference represents more than one order of magnitude in the probability densities. 
The coupling between the two LFD amplitudes has a very small (0.1%) and attractive effect in 
the binding energy. 

In the strong binding limit (B=0.5) the situation is quite similar with enhanced relativistic 
effects. One has azj?D=2.44 for a^R—l.ll and the differences in the wavefunctions - displayed 
in figures [18] and [j~9]- are already visible at k = (figure |18|). One can see however in figure [19] 
that - even for deeply bound systems - fx component still dominates over f 2 . To complete these 
results, we have displayed in figure [17| the LFD, NR and WC coupling constants for different 
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values of the binding energy. One can see that the LFD results are systematically closer to the 
non relativistic values than awe are, as if the fermionic character of the constituents generates 
closer binding energies to the NR case, but larger differences in the high momentum components 
of the wavefunction, due to the different asymptotics of interaction kernels. 




0.5 1 1.5 2 2.5 3 1 2 3 4 5 6 



k k 



Figure 18: LFD and non relativistic wave- 
functions for J=0 + state with B = 0.5, 
fi = 0.15 in Yukawa model 



Figure 19: LFD and non relativistic wave- 
functions for J=0 + state with B = 0.5, 
\x = 0.15 in Yukawa model 



Pseudoscalar coupling 

For pseudoscalar coupling the situation is more involved. First the use of vertex form 
factor - though not necessary for getting convergent solutions - is essential to ensure physically 
meaningful results. We have chosen, as in [0 , a form factor 

with A = 1.3 and n = 1. 

In the weak binding limit (B=0.001) one has cxlfd = 190 and oinr = 166, a repulsif effect 
much stronger (15%) than in the scalar coupling. Corresponding wavefunctions are shown 
in figure One can see that f% ~ /i at k ~ 0.3 and dominates above k=l. It is worth 
noticing the dramatic influence of the form factor in all these calculations: one has for instance 
«Li?D=103 for A = 5 and a^FD=Y12h for A=0.3 ! 

The coupling between the two components is also very important. By switching off 
the non diagonal kernels Byi = -B21 = the coupling constant moves from oilfd = 190 to 
Q-lfd = 251. It has thus an attractive effect which tends to minimize the difference between 
LFD and NR results. 

Quite surprisingly, in the strong binding limit (B=0.5) we have found «lf£i=1462 and 
aivfl=3065. Relativistic effects become now strongly attractive ((Xlfd < &nr)- An essential 
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part of this attraction is due to coupling fx—fa of the two components in the LFD wavefunction. 
By performing one channel calculations one has indeed ajjrD = 3001, what represents a strong 
reduction in the effect though it remains slightly attractive. We have checked if this happens 
for different values of the exchange mass \i. For the same binging energy B = 0.5 and \i — 0.5 
one has aLFo=1728 and 0^^=1400, repulsive once again. This told us the difficulty of talking 
about relativistic effects in general. They turn to depend not only on the kind of coupling 
but also on the binding energy of the system and, furthermore, on the mass of the exchanged 
particle. 

The preceding results show a qualitative difference between S and PS cases. Pseudoscalar 
coupling is by construction relativistic: small and large spinor components are mixed to the 
first order. Moreover the coupling between fx and fa is essential even for very weakly bound 
systems. It is so interesting to study in this case the zero energy limit of the LFD and compare 
with the non relativistic results. We understand by that the static potential, with the same 
form factor, including delta function term. The results, displayed in figure |2T|, correspond 
to \i — 0.5 and two different cutoff parameters A in the form factors. They show the same 
behavior that was found in the scalar case, i.e. that relativistic and non relativistic approaches 
did not coincide even for systems with zero binding energies as fas as they interact with massive 
exchanges. 
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Figure 20: LFD and non relativistic wave- 
functions for J=0 + state with B = 0.001, 
H = 0.15 and PS coupling 



Figure 21: Zero binding limit in LFD (solid) 
and Schrodinger (dashed) equations with PS 
coupling and different form factors (fi = 0.5) 



5 Summary 

We have presented the main ideas of the Light-Front Dynamics formalism in its explicitly covari- 
ant version. When applied to nucleon-nucleon system with perturbative wavef unctions, they 
led to results which are in good agreement with the TJNAF data, even for momentum trans- 
fers larger than the constituent masses. Recently, an interesting application of the Light-Front 
nucleon-nucleon wavefunctions was done to successfully describe the two-body correlations in 
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light nuclei |8j . These facts indicate the ability of this approach in describing composite rela- 
tivistic systems and claim for further developments on this theory. 

This work was pursued by calculating the Light-Front solutions in the ladder approxima- 
tion for two-scalar and two-fermion systems and the main results obtained up to now have been 
summarized in this contribution. Among them we would like to stress the points that follow. 

In the scalar case, we have found that the inclusion of relativity has a dramatic repulsive 
effect on binding energies even for systems with very small values. The effect is specially 
relevant when using a scalar model for deuteron: its binding energy is shifted from 2.23 MeV 
down to 0.96 MeV. This can be corrected by a 10% decrease of the repulsive coupling constant, 
what indicates the difficulty to determine beyond this accuracy the value of strong coupling 
constants within a non relativistic framework. 

Light-Front wavefunctions strongly differ from their non relativistic counterparts if they 
are calculated using the same coupling constant. Once the interaction parameters are readjusted 
to get the same binding energy, both solutions become closer in the small momentum region 
but their deviations are sizeable at k ~ m. 

The relativistic effets are shown to be induced mainly by the additional terms appearing 
in the interaction kernel. Kinematical corrections have only a small influence on the binding 
energy. 

The Light-Front results are found to be quite close to those provided by Bethe-Salpeter 
equation, for a wide range of coupling constant, despite the different physical input in their 
ladder kernel. 

Even in the zero binding limit, Light-Front and Schrodinger solutions differ (/j, ^ 0) both 
for scalars and fermions. This leads to the conclusion that such systems cannot be properly 
described by using a non relativistic dynamics. 

For the Yukawa model, we found a critical coupling constant below which the J=0 + 
solutions are stable without any regularization. For the pseudoscalar coupling, these solutions 
are also stable for all values of the coupling constant but lead to a quasidegenerate a(B) relation. 
The J=l + states are on the contrary unstable for both couplings. 

Pseudoscalar coupling shows large deviations with respect to non relativistic solution. 
Wavefunction is dominated by relativistic components at k « m, even for weakly bound 
systems. However form factors play there a determinant role. 

The question about relativistic effects has no simple answer. The consequences of imple- 
menting the Lorentz invariance in a quantum mechanical description of a system are different, 
following: the nature of its constituents, the kind of interaction, the quantum numbers of the 
state, its binding energy, and even the mass of the exchanged particle! There are no simple 
recipes to perform a priori evaluations. 
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